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Sound  Propagation  in  the  Sea  — Ray  Tracing 

CURVATURE  CORRECTIONS  FOR  THE  EARTH  (SEA  LEVEL)  SHAPED  AS  AN 
ELLIPSOID  OF  REVOLUTION  k*  i 

by:  M.  M.  Holl 
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The  Problem 


The  three-dimensional  tracing  of  a sound  ray  in  the  sea  may  be 
referred  to  the  surface  coordinate  system  of  latitude  and  longitude  and  to 
the  depth  below  the  sea  surface.  This  coordinate  system  involves  curva- 
ture terms . t’1' v 

b*  (,i 

■ For  the  earth  shaped  as  an  ellipsoid  of  revolution  may  refer' to 
two  definitions  of  latitude:  the  geocentric  latitude  <t>  and  the  astronomical 
latitude  $ . These  are  illustrated  in  Fig.  2.  There  is  a one-to-one 
relationship  between  these  definitions  and  the  transformation  is  simple  as 
will  be  shown.  W*-eh©o&»  to  e*<*  the  geocentric  latitude  $>  ips  principal 
latitude  coordinate.  7*' 

Let  the  direction  of  propagation  of  a sound  ray  at  an  arbitrary 
point  in  its  progress  --  latitude  <p  , longitude  and  depth  z --  be  defined 
by  the  unit  vector,  T . It  makes  the  angle  & with  the  horizontal  as  shown 
in  Fig.  1.  The  horizontal  component  of  the  unit  vector  T defines  the 
direction  of  the  unit  vector  It  . 

Vi  v\ 
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Fig.  i Definition  of  Parameters  In  the  Vertical  Plane  of  the  Kay 

In  Fig.  2 we  define  and  illustrate  geometric  parameters  which  are 
relevant  to  the  problem.  These  Include  the  angle  jA  which  It  makes  with 
the  unit  vectoi  i : t is  directed  £ radians  north  of  due  east. 

The  coordinate-curvature  terms  arise  In  prescribing  the  ray 
direction  IT  by  the  angles  ^ and  & . I he  angle  ^ is  ie toned  to  the 
horizontal  which  Itself  is  turning  in  space  from  location  to  location.  The 
angle  S is  referred  to  the  east  which  itself  is  turning  from  location  to 
location. 

We  may  formulate  the  curvature  terms  directly,  lhe  reference 
directions  turn  only  by  horizontal  progress  and  the  angles  are  defined 
in  specific  planes.  The  terms  are 
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Fig.  2 Definition  of  Parameters  for  the  Surface 


.-i 


$ : geocentric  latitude 
^ : astronomical  latitude 

0<  « $ 


1^  : unit  vector,  surface  normal , 
alligned  with  local  gravity, 
subtends  angle  with 
equatorial  plane  as  shown. 

4 : unit  vector,  surface  tangent, 

directed  eastward. 

J : unit  vector,  surface  tangent, 
directed  northward. 

i x j > Ik 


A 

B 

E 

m , h 


t 

in 


equatorial  radius  of  earth  A > B 
polar  radius  of  earth 

radius  at  geocentric  latitude  <p  (astronomical  latit  de  ) 

Cartesian  coordinates  in  arbitrary  meridional  plane  of  longitude  9 . 

Solid  lines  lie  in  the  meridional  plane;  dashed  lines  do  not. 

unit  vector,  tangent  to  surface  trace, 
unit  vector,  surface  tangent,  normal  to  it  . 

Ik  x ft  * in 

t = cos  /8  i + sin  $ 

tn  - - sin  A a + cos  ^ 


N 


| 

! 
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2. 


Geometry  of  the  Surface 


The  intersect  of  the  arbitrary  meridional  plane  and  the  surface  is 
an  ellipse: 


+ 


2S 


1 


(3) 


We  substitute 


h - E sin  $ 
m = E cos  <$> 


(4) 


and  obtain 


A 


B 


2 2 2 2 , 
A sin  c|>  + B cos  <£> 


(5) 


We  also  note  that 


tan  <p 


h 

m 


(6) 


An  incremental  northward  displacement  on  the  ellipse  may  be 
expressed  by 


h^  ♦ m**  - 0 (7) 

A 
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where  &h  is  positive  and  Sm  is  negative.  This  displacement  is 
directed  along  | . Hence 


tan  $> 


% m _ A^  h 

S h B2  m 


A 


tan  <|> 


(8) 


This  is  the  relationship  between  the  astronomical  latitude  $ and  the 
geocentric  latitude  <$>  . Their  difference  angle 


oc  * ^ - 4> 


(9) 


is  given  by 


tan  «■  * tan  ($-$>) 


tan  $ - tan  » 
1 + tan  $ tan  4> 


(A^  - B^)  tan  <fr 

2 2 2 
B + A tan  $ 


We  also  require  the  ratio  of  latitude  increments:  Differentiation 
of  Eq.  (8)  yields 


2 < c -f  A 2 . c . 

sec  ^ = — r sec  4>  dip 

B 


_S|_ 

S $ 


2 i 

cos  $ 

2 * 
cos  q> 


A2  B2 


2 . 4 4 2 v 

cos  4>  (B  + A tan  <J>) 
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3.  Curvatures  of  the  Coordinate  Unit  Vectors 

We  are  now  in  position  to  express  and  develop  the  curvatures  of 
the  coordinate  unit  vectors.  These  are  expressed  by  an  incremental 
directed  turning  of  the  unit  vector  divided  by  the  space  increment,  in  the 
pertinent  direction,  over  which  this  turning  takes  place: 


l ■ v i 


E cos  <p  JM' 


( i sin  £ - Ik  cos  $ ) 


(12) 


I 


i sin  $ Jte 

E cos  4 


(13) 


i cost 
E cos4> 


(14) 


i • v £ 


0 


i ’ 


k 6$ 

E sec  « 


(15) 


(16) 


4 ' ^ ^ E sec  o,. 


Ik  • V 1 


Ik  • v4  = Ik-  v K = 0 


In  Eqs . (16)  and  (17)  we  substitute  from  Eq.  (11): 


i -vi 


. .2  n2 

R COS  e(  A B 

E cos 2 4>  (B4  + A4  tan24>  ) 


. v fc 


A A2  o2 

a cos  tx  A B . 

' 7 4 4 2 v 

E cos  (B  + A tan  <J>  ) 


npyplripment  of  the  Curvature  Terms 


Equation  (1)  is  developed  as  follows 


( d i / b s) 


t • V t • fe 


(cos /9  * + Sin£  £)  • v (cos£  a + sin(J  £ ) ■ k 


“g  g . v ( • Ik  + cos  fi>  sin/3  « • ^ -J 


j . v • Ik 


+ sin  jft  cos  /3  i • v i - k + sin2/3  | • * { . k 


cos  cos 
E cos  4> 


2 2 2 
sin  4 cosot  A B . 

E cos24>  (B4  + A4  tan2  4>  ) 


Equation  (2)  is  developed  as  follows 
(i/9  /a  s)c  = - t • v i ■ { 


- (cos  4 i + sin  £ ^ ^ i • -J 

- cos/J  i v i .j  - sin  £ j-vi  • 


cos  >6  sin  $ 
E cos  $ 


(2  2) 


We  now  wish  to  eliminate  E,  f and  * by  substitution  in  terms  of 
$ according  to  Eqs.  (5)  , (8)  and  (10).  For  abbreviation,  where  convenient, 

we  adapt 


K = tan  <$>  ; F = 9 

B 

We  may  rewrite  Eqs.  (5),  (8)  and  (10)  as  follows: 

2 = A2  (1  t lA 

(1  + FK2) 

tan  § = FK 

(F  - 1)  K 
ts  n “ o 

(1  + FK  ) 


(23) 


(24) 


(25) 

(26) 
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It  follows  that 


cos  $ = 


2 2 

(i  + r o 


sin  £ * 


cos  c*.  = 


2 2 

(i  + r n 


(1  + F K ) cos 

2 2 1/2 

a + f2  n 


With  the  appropriate  substitutions,  Eqs.  (21)  and  (22)  become 


(i«/is) 


( d/3/a  s) 


1 1 1 + fk 

A ) , . 


1 + F2  K2 


cos  A 
A 


cos2#  + r sin2l  + r2K2 

2 2 
1 + F K 


1 + F K 
1 + F2  K2 


These  are  the  coordinate-curvature  terms , on  and  (&  , as  functions  of 
the  geocentric  latitude  ^ anc^  ^he  ray-direction  angle  /3  . 

In  the  special  limit  of  a sphere: 

A * B = R;  4>  = § ; F = 1 


the  curvature  terms  simplify  to 
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( wa  s) 


i 

R 


(32) 


C,R 


( i Z3  / i s)c  R 


cos  d tan  4> 
R 


(33) 


5 . The  Ray  Tracing  Equations 

The  ray  tracing  equation,  in  vector  form,  may  bo  expressed  by 

T • v T - T x ( T x v m C)  (3-1) 

where  T is  the  unit-vector  ray  direction  and  C is  the  sound  speed.  We 
take  the  components  along  N and  m — unit  vectors  defined  in  Tigs.  (1) 
and  (2)  respectively  — and  obtain 


T • 

v T • N * 

- IN  • V In  C 

(35) 

T • 

V T • m 

- H • V InC 

(36) 

The  left-hand  sides  may  be  transformed  by  geometrical  consideration: 


T 

• V T . 1 N * 

i J / J S - cos  if  t • v t • Ik 

(37) 

T 

• v T • n = 

cos  yj  Was  + cosy  t - v i - j] 

(38) 

where  S is  a linear  measure  along  the  ray  trace: 

§ s = S S cos  V (39) 


- 
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The  coordinate-curvature  terms  are  expressed  by  Eqs  . (1)  and  (2) 
and  the  ray-tracing  equations  may  be  written 


sec  ^ IN  • V In  C 


# \n  • V In  C 


The  curvature  terms  are  given  by  Eqs.  (30)  and  (31) . 


Linearization  as  to  Oblateness  of  Earth 


If  we  define  € by 


(1  + € ) B 


where  A is  the  equatorial  radius  and  B is  the  polar  radius  of  the  eaith, 


0 < £ < < 1 , H 

and  we  may  linearize  Eqs.  (30)  and  (31)  by  dismissing  squares  and  higher 
powers  of  £ as  negligible.  In  this  linearization 

A2  = (1  + 2€  ) B“  (4 


(1  + 2 € ) 


(1  + 4€  ) 


We  also  define  as  the  earth's  mean  radius 


It  follows  that 


R = 2 (A  + B) 


(1  + 2 € } R 


With  these  substitutions  we  linearize  Eqs.  (30)  and  (31)  and 


obtain 


(ft)  --if 

c v. 


1 ( 2 2 2 1 
— ( 1 + € (2  sin  /3  cos  <J>  - sin  <f)  - ^ ) 


Or)  - - 


cos  $ tan  <t> 
R 


1 + 6 + cos2  ) j 


We  note  that  as  € — * 0 we  arrive  at  Eqs.  (32)  and  (33)  respectively,  for 
the  sphere;  this,  however,  is  not  a check  on  the  linear  terms. 

It  is  convenient,  at  this  point,  to  note  that  the  effect  of  the  depth 
of  the  ray  on  the  curvature  terms  can  be  taken  into  account  by  replacing  R 
by  R - z.  This  is  equivalent  to  replacing  A and  B by  A - z and  B - z in 
all  earlier  expressions  of  the  curvature  terms;  the  effect  on  e is  negligible, 
In  general  the  depth  correction  is  smaller  than  the  oblateness  correction. 


-12- 


For  convenience  we  also  exhibit  Eqs.  (49)  and  (50)  in  terms  of  the 
astronomical  latitude , $ : 


By  linearization  of  Eq.  (26)  we  find  that  the  difference  between  the 
astronomical  latitude  <|>  and  the  geocentric  latitude  $ is 

<x  = $ - <£>  = 6 sin  2 

We  note  that 

oc  = sin  <x  =»  tan  « = fe  sin  2 <f>  = € sin  2 $ (53) 

cos  <x  = sec  <x  = 1 (54) 

The  angle  tx  Is  a maximum  at  about  <$>  * + 45°. 


Magnitude  of  Corrections  for  Oblateness 


The  oblateness  factor  does  not  enter  in  sea-surface  or  sea-bottom 
reflections.  These  reflections  are  treated  in  the  same  way  as  in  the 
spherical  case;  no  oblateness  correction  enters;  this  is  because  the  angle 
t is  defined  relative  to  the  local  vertical  (see  Fig.  1). 
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The  values 


g 

A - 6.378  x 10  meters 
B * 6.357  x 10  meters 

have  been  quoted;  perhaps  better  values  are  now  available.  These 
give 


r - 1 « 0.0033 


300 


The  oblateness  corrections  have  the  following  ranges  in  the 
ture  terms: 


(ft),  --if 


1 + (-  0.005  to  + 0 


.005)j 


for  both,  Eqs.  (49)  and  (51). 


(ft) 


cos  tan 
R 


— | 1 + (0.0017  to  0.005)] 


| ilj  * _ cos  <3  tan  j l + (_0>005  to  -0.0017)] 


for  Eqs.  (50)  and  (52),  respectively. 

In  moving  northward  at  the  equator,  we  have 


( **/*  s) 


- ( 1 + 0.005) 

K 


(55) 

values 


(56) 


curva- 

7 

(57) 


(58) 


(59) 


(60) 


A j 
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In  a northward  displacement  of 


b s = 0.1  R = 637,000  meters 


(61) 


the  oblateness  correction  amounts  to  0.0005  radians.  Spreading  this 
correction  linearly  over  the  displacement  results  in  a depth  correction  of 


637,000  x 0.00025  = 160  meters  (62) 

This  is  an  upper-limit  approximation  to  the  displacement  errors  made  in 
neglecting  the  earth's  oblateness. 

We  also  note  from  Eg.  (53)  that  maximum  c*  (at  ~ i 45°)  is 


0.0033  radians  ^ 11  minutes 


(63) 


8 . The  Geodesic:  Surface  Ray  with  no  Refraction 

Consider  a ray  which  is  trapped  in  the  sea  surface  by  continuous 
glancing  reflection,  but  undergoing  no  refraction.  Its  trace  describes  a 
geodesic: 

ft  . v t - (t  . v ft  • Ik  IV*  (64) 

Its  geodesic  curvature  (l.e.  curvature  in  the  surface)  is  zero: 

t • V t . in  * 0 (65) 

We  have  seen  that  Eg.  (65)  transforms  into 

bfi/bs  - (i*/*s)  (66) 


- L* 
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We  choose  Eq.  (22)  to  write 


t 


^ B _ cos  0 sia  $ 

is”  E cos  <J> 


for  the  governing  equation  for 


a 


geodesic. 


(67) 


The  governing  equation  can  be  integrated  once  and  obtains  that, 
for  a geodesic,  /3  is  a function  of  latitude  only.  We  achieve  this  by 
substituting 


sin  ft  S s = E sec  <x  & 4> 


(68) 


as 


seen  from  Fig.  2,  into  Eq.  (67)  to  eliminate  & s.  We  obtain 


tan  ft  S ft  = 


in  $ 


sin 


cos  ()>  cos  <*. 


S 4> 


(69) 


We  transform  Eq.  (28)  to  obtain 


V 


sin  $ = sin  4> 


| 1 + (F2  - 1)  sin2<f>] 


1/2 


(70) 


We  transform  Eq.  (29)  to  obtain 


cos  «. 


1 + (F  - 1)  sin  <t> 


| 1 + (F2  - 1)  sin2 


* 1 


1/2 


Introduction  of  Eqs.  (70)  and  (71)  in  (69)  yields 


tan  ft  %ft  = - tan  <}> 


1 + (F  - 1)  sin  <$> 


S <$> 


(71) 


(72) 


-16- 


We  integrate  Eq.  (72)  from  the  equator,  4>  ■ 0,  where  /3  - ^o, 

to  an  arbitrary  latitude  $ . The  result  is 


2 1/2 

cos  /3  =>  cos  /3  (1  + T tan  <$>  ) 


The  maximum  latitude  which  this  ray  can  attain  is  reached  when  /3  = 0 : 


tan  4) 


tan^  A, 
F 


The  analytical  equation  for  the  geodesic  requires  one  more  inte- 
gration. We  first  transform  Eq.  (73)  to  obtain 


tan  = 


2 o 
sec  '3o 


1 + F tan  <)> 


From  the  geometry  of  Fig.  2 we  obtain 


tan  /9 


sec  <*  S 4> 

J?  cos  4>  S 9 


It  follows  that,  for  the  geodesic, 


d ♦ . \ sec 

— r = COS  (f>  COS  01  \ “-T— 

d 0 1 1 + F tail  4> 


Equation  (77)  may  also  be  expressed  by 


t 


d 0 


^ 0 1/2 

(2 
sec  <t>  sec  <*  ( 1 F tan"  ft  ) 


J 


1/2 


(tan  B 


T tan*"  ft  ) 


d.  ft 


(78) 


The  geodesic  trace  is  obtained  by  4>  cycling  back,  and  forth  between  ftm 

and  - ft  , the  range  given  by  Eq.  (74). 
r m 

In  the  linearization  of  Eq.  (78),  as  to  oblateness,  we  must  concern 
ourselves  with  the  range  of  ft  . We  may  immediately  set  sec  oc  equal  to 
one;  its  departure  is  second  order  in  € . Again  using  K 5 tan  ft  we 
transform  Eq.  (78)  into 


K 

/ 

2 1/2 

9 ' 9o  * \ 

(1  t 

O 

F K1 

1/2  2 

dK 

2 1/2 

(79) 

J 

0 

(1  + O 

(tan  B 

0 

fk  ) 

We  change  the  variable 

of  integration  once  1 

more 

to 

P 3 

f1/2k 

(80) 

and  obtain 

P 

( 

F~  1/2  fl  + P2) 

1/2 

d P 

(81) 

9 - 9 - \ 

° J 

0 

(1  + 

-1  2 1/2 
r p ) 

- 1/2 

(tan2  3 - P‘) 

0 

Linearization,  following  F ■ 

1 + 2 € , yields 
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p 


p 


9 - 9o  = \ 2 2 1/2 

oJ  (tan2/3o-P2) 


J (1  + P2)  (tan2  /3  - P2) 
o o 


where  the  range  of  P is  _+  tan  /?  . 

We  introduce  to  by 


sin  to  = P cot  /3 


2 2 
tan  to  = P (tan  /3q  - P ) 


Equation  (82)  transforms  into 


arcsin  (P  cot  /3  ) 
, o' 


0-0  = 

o 


arcsin  (P  cot  (3  ) 
o 


2 2 
1 + tan  <3  sin  u> 
o 


The  first  integral  is  solved;  the  second  is  transformed  once  more: 


P j tan2  3 - P2} 

i o J 


9-0  = arcsin  (P  cot  /3  ) - € 

o o 


J 1 + sec2  3 y2 
o o 


where  y replaces  the  dummy  variable,  tan  <o 


' 


The  linearized  solution  for  the  geodisic  is 


0-0  = arcsin  (P  cot  6 ) 

o o 


{ 


2 2 

€ cos  /3  arctan  ^ sec  ft  P (tan  ft  - P ) 
o I o o 


•1/2  j 


(86) 


where 


P = (1  + t ) tan  $ 


(87) 


has  the  range  + tan  Q . The  longitude  at  the  maximum  latitude  is 


0 - 0 = V2  (1  ~ « cos  fi  ) 

mo  o 


(88) 


For  £ = 0 the  surface  becomes  a sphere,  and  the  geodesic  is  a 
great  circle:  Equation  (86)  reduces  to 

sin  (0  - 0q)  = tan  4>  cot  j3Q  (89) 

This  result  could  have  been  obtained  directly  from  the  spherical  trigo- 
nometry of  a right  spherical  triangle. 

In  Fig.  3 we  show  the  extent  to  which  orbits  of  the  geodisic 
precess.  For 


£ 


1 

300 


and 


/3  = 45  degrees 


(90) 
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Fig.  3 Schematic  ot  the  Geodesic 


the  latitude  maximum  attained  in  Fig.  3 is  located  at 


0 

89. 7g 

degrees 

4>  - 

44. 90 

degrees 

(^)  5K 

45. 10 

degrees 

(91) 

The  oblateness  correction  appears  to  be  negligible  within  the  practical 
limitations  in  sound-ray  tracing  in  the  sea. 


